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ABSTRACT 

Models for the characteristically soft X-ray spectrum of the compact X-ray source V1408 Aql 
(=4U 1957+11) have ranged from optically thick Comptonization to multicolor accretion disk models. 
We critically examine the X-ray spectrum of V1408 Aql via archival Advanced Satellite for Cosmology 
and Astrophysics (ASCA) data, archival Rontgensatellit (ROSAT) data, and recent Rossi X-Ray Timing 
Explorer (RXTE) data. Although we are able to fit a variety of X-ray spectral models to these data, we 
favor an interpretation of the X-ray spectrum as being due to an accretion disk viewed at large inclination 
angles. Evidence for this hypothesis includes long term (117 day, 235 day, 352 day) periodicities seen by 
the RXTE All Sky Monitor (ASM) , which we interpret as being due to a warped precessing disk, and a 
1 keV feature in the ASCA data, which we interpret as being the blend of L fluorescence features from 
a disk atmosphere or wind. We also present timing analysis of the RXTE data and find upper limits of 
4% for the root mean square (rms) variability between f=10 — 3 — 16 Hz. The situation of whether the 
compact object is a black hole or neutron star is still ambiguous; however, it now seems more likely that 
an X-ray emitting, warped accretion disk is an important component of this system. 

Subject headings: accretion, accretion disks — black hole physics — neutron star physics — stars: 
individual (V1408 Aql) 



1. INTRODUCTION 



clination of i ~ 70° -75°, consistent with the models for 



The low mass X-ray binary (LMXB) V1408 Aql 
(=4111957+11, 3U 1956+11) was detecte d during scans 
of th e Aquila region by Uhuru in 1973 flGiacconi ct aT 



1974 ), and it was subsequently i dentified with an 18'.° 7 
star having a s trong blue excess ( Margon, Thornstcnscpj 



& Bowyer 1978). The object is situated in a region ofrel- 

; 1.3 x 10 21 
1992|) 



atively small exti n ction (Yh 
& Lockman 199C; stark et al 



cm , Dickey 



Ay measurements 
place the source at a distance >2.5kpc, and comparisons 
of its X-ray and optica l luminosity to Sco X-1 place it at 
a distance of ^7kpc (Margon, Thornstensen & Bowyer 



1978) 



Little is known about the nature of the system. Optical 
spectra of V1408 Aql reveal a power-law continuum with 
Ha, H/3, and Hen 4686A emission lines (Cowley, Hutch-| 
ings fc Crampton 1988 : Shahbaz et al. 1996[), typical for 



the optical variability (Hakala, Muhli & Dubus 1999). 

V1408 Aql is one of the less well-studied possible black 
hole candidates (BHCs). Identification as either a BHC 
or a neutron star-low mass X-ray binary (NS-LMXB) is 
usually made by analogy with the spectral- and timing- 
behavior of better observed sources. V1408 Aql has been 
a BHC since 1984, when EXOSAT X-ray observations re- 
vealed that V1408 Aql has a very soft X-ray spectrum, 
similar to that of other BHC. In color-color-diagrams, 
V1408 Aql lies halfway between the black hole candidate 
GX 339—4 (in its high/soft stat e) and the neutron-star 
LMXBs Cyg X-2 and LMC X-2 flWhite fc Marshall 1984 ; 
Schulz, Hasinger fc Triimpcr 1989| ). This color identifica- 



an accretion disk-dominated system. |Lhorstcnscn (1987)| 



reported a nearly perfectly sinusoidal V - band luminosity 
modulation with 10% amplitude and a 0.380 d (— 9.33h) 



tion of V1408 Aql as a BHC, however, is not definitive 

Previous narrow-band observations have not character- 
ized the X-ray spectrum in a consistent manner. The anal- 
ysis of 1983 and 1985 EXOSAT observ ations of V1408 Aql 
led to contr adictory results. While |Singh, Apparao ~& 



Kraft (1994) succeeded in fitting a Comptonization spec- 



period, which he interpreted as due to X-ray heating of 



the companion. — In recent multicolor photometry a more 



trum to these data and interpreted this a s an indication 
that V1408 A ql is a black hole candidate, Ricci, Israel & 



, . , , ' .,, 0(W j , , . ,., j Stella (1995) interpreted the same data as being similar to 

complex hghtcurvo with 30% modulation amplitude wao -A—— A T A/rY „ n , ,„„ * , n . 



observed. Hakala, Muhli fc Dubus (1999)| interpret this 
change in the shape of the lightcurve as evidence for a disk 



with a large outer rim, possibly due to a warped d isk, accn 



elope to edge on (Hakala, Muhli fc Dub up 1900; ooo aloo 



that observed from NS-LMXBs. Observations with Ginga 
with its larger spectral range and effecti ve area, have shed 
more light on the nature of V1408 Aql ( Yaqoob, Ebisawa 
fc Mitsuda 1993). The values of the normalizations of 



. This interpretation is also consistent with the sh ape of 
the infrared spectrum (Smith, Beall & Swain 1990| ). The 
short orbital period is indicative of a late type main se- 
quence star of M rv 1 M m as the donor star — The absence 



multicolor disk models (MCD; [Mitsuda et al. 1984IJ , i.e. 
(r- ln /d) 2 cos i where r ln is the inner disk radius, d is the dis- 
tance to the source, and i is the inclination, have been used 



Of X-ray P plip S P S and thp agsiimpt.in7 t.Tmt. thfl donor star Lcwin 1995Q . In the case of V1408 Aql 



to distinguish between BHCs and NS-LMXBs ( |Tanaka & 

cos 1 / 2 i ~ 2 km 



fills its Roche lobe yield an upper limit on the orbital in- 



assuming d — 7kpc, which is more characteristic of sources 
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V1408 Aql (4U1957+11) 



containing neutron stars. Additionally, the Ginga observa- 
tion showed evidence of a hard tail (1-18 keV) comprising 
~ 25% of the inferred flux for this system at that time. 
The best fit power-law photon indices for the hard com- 
ponent ranged from T ss 2 to 3. 



The EXOSAT observations of |Ricci, Israel fc Stella 
(19950 indicate the presence of an iron fluorescence line 
with an equivalent width of 90 eV or smaller and a line- 
energy of 7.06 keV (i.e., highly ionized). Ot her values in 
the li terature range fr om non-det ection (e.g., Yaqoob, Ebi 



sawa Mitsuda 1991 ) to 200 eV flWhite fe Marshall 1984^ 



the uncertainty being mainly due to the difference in the 
assumed spectral continua and the different sensitivities of 
the instruments. 

Except for one observation, which hints toward a weak 
red-noise (/ _Q ) component between 10~ 4 and 10 _3 Hz, all 
EXOSAT observations are consistent with the abs ence of 
any periodic features (Ricci, Israel & Stella 1995). The 



Ginga observations have yet to have their short timescale 
variability analyzed; however, they do show evidence of 
significant flux and color changes on long time scales 
(> 10 4 s). The Vela 5B satellite did not detect any long- 



term X-ray variability from this source ( Priedhorsky & 
Terrel p984 ) ; however, the upper limits to the variability 



were $loi particularly strong. 

tie JJUblMlUd ypemi'al lllOdely are 1 aWieTjle'd al faCe 1 
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2. SPECTRAL ANALYSIS 

V1408 Aql was observed with RXTE in 1997 Novem- 
ber in three observing blocks for a total on source time 
of 27ksec. A log of the observations is given in Table 0. 
Since the spectral shapes of the three observing blocks are 
identical, the data were analyzed together. Spectral and 
temporal dat a we re extracted using the methods outlined 
in appendix A.l. Spe ctral analysis was performed with 
XSPEC, version 10.0s ( [Arnaud 199^ ). 

The RXTE spectrum of V1408 Aql is very soft. The 
Proportional Counter Array (PCA) did not detect any flux 
above ^20keV and the High Energy X-ray Timing Exper- 
iment (HEXTE) count rates are consistent with zero: the 
background subtracted count rates were 0.4 ± 0.2 cps and 
0.0 ± 0.2 cps, for HEXTE clusters A and B, respectively. 
The residual flux in cluster A is most probably due to 
a slight overestimation of the HEXTE background dead 
time, as the spectrum seen is similar to the HEXTE back- 
ground. Therefore, we do not consider V1408 Aql to be 
detected with HEXTE and will not further discuss these 
data. 

To describe the PCA spectrum we use the spectral mod- 
els traditionally applied to V1408 Aql: an expo nentially 
cutoff pow er-law, a multicolor disk-black body (jVIitsud 
et al. 1984), and a Comptonization model after Titarchu' 



(1994) . Due to the low sensitivity of the P CA to the low 



value 1 bfaErr-tfae y dulive energetics of the- disk black-body absor bing column towards V1408 Aql (see |Stelzer et al 



and 



nower-iaw components, as well as xhe slope of the 
high-energy power-law, are very similar Lu Lhoso seen in 



BHCs such as LMC X-l ( Ebisawa, Mitsuda fc Inouc 1989 
Wilms et al 1998b| ), LMC X-3 ( fTreves et al. WM iWilm' 



998b[ ) and in the soft state of GX 339-4 (|Miyamotc 
1991; prcbenev et al. 1991). However, at luminosi- 



et al. 
et al. 

ties as low as that of V1408 Aql, BHCs tend to show hard 
tails with no evidence of a disk or thermal component. On 
the other hand, NS-LMXBs that exhibit soft disk spectra 
also tend to show an additional ~ 2keV blackbody com- 
poncnt, while showing little hard flux ( Miyamoto 1994 , 
and references therein). Furthermore, low-luminosity NS- 
LMXBs also tend to be dominated by hard emission. 

Thus, there are good arguments that point towards 
V1408 Aql being a neutron star and also toward it being a 
black hole; however, none of the arguments are truly con- 
clusive. In either case, V1408 Aql would still be a unique 
object, being either an unusually soft low-luminosity BHC, 
an unusually soft low-luminosity neutron star, or a soft 
neutron star with an unusually energetic hard tail. With 
the advent of X-ray detectors with much larger effective ar- 
eas than EXOSAT and Ginga, as well as with the availabil- 
ity of detectors of higher energy resolution, such as those 
on the Advanced Satellite for Cosmology and Astrophysics 
(ASCA), a critical reexamination of the X-ray spectrum of 
V1408 Aql has become possible. In this paper we present 
the results from our analysis of a 30 ksec pointed observa- 
tion with the Rossi X-Ray Timing Explorer (RXTE), as 
well as archival data from ASCA and the Rontgensatellit 
(ROSAT). In §|| we present the results from the spectral 
analysis. We discuss the timing analysis in §[|, the long 
term variability of the source in §|], and we discuss our re- 
sults in §||. The details of the data extraction are described 
in an appendix. 



1999 for a discussion of the sensitivity of the PC A to 
Ah), we fixed Ah to the pickey fc Lockman (1990") value 
of Ah = 1.3 x 10 21 cm . The results of our spectral 
fits are given in Table [| while the PCA spectrum and 
the residues are displayed in Fig. [|. All three models 
roughly describe the observational data. Note that we 
do not see any evidence for a high energy power-law tail 
as that seen in previous observations. The 90% confidence 
level upper limit to the 3-20 keV flux from a power law 
is 8 x 10 -12 ergs cm' 2 s _1 , which is less than 2% of the 
observed 3-20 keV flux. 

The best description of the PCA data is given by the 
Comptonization model (x 2 c d = 0-82 for 29 degrees of free- 
dom), while the residues of the MCD model and the ex- 
ponentially cut-off power-law show structure in excess of 
that expected from calibration uncertainties of the PCA. 
These residues are especially apparent in the low energy 
channels of the PCA, below the characteristic feature of 
the Xe L-edge at ~ 5 keV (a region of very uncertain detec - 
tor calibration; see the discussion by Wilms et al. 1998a). 



Inspection of our best fit values in Table |2j shows that the 
Comptonization model results in such a good fit because 
the seed photon temperature of the model, taken here as a 
Wien spectrum with best-fit temperature kTo = 0.34 keV, 
is uncharacteristically large. This is also evident in the 
very asymmetric confidence contour indicated in Table ^. 
Setting kT to a small value yields residues that more re- 
semble those seen in the MCD and exponentially cut-off 
models. We conclude that there is unambiguous evidence 
for the presence of an additional soft-excess below ~ 5 keV. 
Since the low-energy cut-off of the PCA is at ~ 2 keV, this 
instrument cannot be used to further constrain the na- 
ture of this soft-excess. We therefore turned to archival 
data from observations of V1408 Aql made with ASCA 
and ROSAT. 
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Table 1 

Observing log for V1408 Aql. Count rates given are background subtracted. 



Date 


On-Sourcc Time 


Instrument 


Count Rate 




[ksec] 




[cps] 


1997 November 26.02-26.45 


21.2 


RXTE-PCA 


230 


1997 November 27.08-27.17 


4.7 




229 


1997 November 29.08-29.15 


1.5 




233 


1994 October 31.13-31.93 


21.6 


ASCA SISO 


9.8 






ASCA SIS1 


7.0 






ASCA GIS2 


19.4 






ASCA GIS3 


21.5 


1992 May 08.39-11.14 


14.6 


ROSAT PSPC 


19.2 




Fig. 1. — Spectral modeling of the PCA data. Residues are shown as the contribution to x- a ) Count rate spectrum and the best 
fit MCD model with an optically thin plasma spectrum, b) Residues of the best fit exponentially cutoff power-law, c) residues of 
the best fit MCD model, d) residues of the best fit Comptonization model, and e) residues of the best fit MCD plus optically thin 
plasma spectrum. 
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Table 2 
Results of spectral fitting 



Model 


Detector 




Parameters 




xVdof 


cutoffpl 

b 

power 


PCA 
SIS 

GIS 
PSPC 


T = -0.26 ±0.06 
T = 0.43 ±0.04 
A/ H = 1.4 ±0.01 
T = 0.56 ±0.02 

T = 1.54 ±0.05 


A PL = 0.21 ± 0.01 
A PL = 0.200 ± 0.002 

yl P L = 0.267 ± 0.001 
A PL = 0.148 ±0.004 


E lo i d = 1.44 ±0.02 
Sfoid = 2.30 ± 0.07 

E loU = 2.44 ± 0.03 
JV H = 1.8 ±0.1 


60.2/30 
559/501 

1630/1342 
56.7/21 


diskbb 


PCA 
SIS 
GIS 
PSPC 


fcTi n = 1.319 ± 0.003 
fcTn = 1.530 ± 0.001 
fcT„ = 1.487 ± 0.004 

fcTin = 1-0 


A«.*m, = 15 7 ± 26 

^ A diskbb J-t/. i — — J ■ 
j4diskbb = 8.7 ± 0.2 

^diskbb = 12.2 ±0.1 

^diskbb = 30.3 


JV H = 1.32 ± 0.04 
A^h = 1.0 


42.0/31 
628/502 
2148/1343 
104/22 


compTT 


PCA 
SIS 

PTC 

PSPC 


fcT = 0.338^o °3 7 
norm = 0.33^;^ 6 
kT = 0.32 ± 0.01 
norm = 0.239 ± 0.003 

l^rji n 07 ±(1(11 

K±o — U.Zl rt U.U1 
norm = 0.333 ± 0.005 
fcT = 0.176 ± 0.030 
norm = 0.15t u ;j£ 


fcT e = 1.17 ±0.01 

kT c 1.28_q Q2 
N H = 0.37±g;JS 

/CJ e — 1.Z4 zt U.U1 

fcT c = 2.40lfo 
N H = 1.07 ±0.30 


r = 10.7 ±0.3 
T = 12.2+°;? 

_ ioo±flO 

T = 9.8i»;| 


23.8/29 
580/500 

lo / 1/ lo41 

29.5/29 


cutoffpl+ray 
power+ray 


PCA 
SIS 
GIS 
PSPC 


T = -Q.69Lto:?g 
fcT ray = 1.17 ±0.19 
T = 0.36 ±0.05 
uqi — i r)8 +0 06 

-L \J.iJfJ —1— KJ.KJO 

A/Tray = 0-^^—0.09 

r = 1.2 ±0.1 

urn _ i nq+ 006 
■Wray — r.u»_ og 


Apl — U.12_ 02 
4 _ 44+23 

A PL = 0.196 ±0.003 
4 _ fi 7 +i.7 

^ray — U. 1 _i,g 

A-dt — 951 ± 002 

/t _ QO + 9 
^ray — tJ ' £ '_3 

4 PL = 0.12 ±0.01 
4 — 1 9 fi+ 4 - 8 

-riray — — 5.2 


£? Md = 1.35 ± 0.06 

Sfoid = 2.21 ± 0.07 
JV H = 1.4 ±0.01 

Ah = 1.3 ± 0.02 


17.4/28 
507/499 
1 303/1 340 
19.1/20 


diskbb+ray 


PCA 
SIS 
GIS 
PSPC 


fcTin = 1.33 ±0.01 
-<4-ray lOO^^g 

fcT„ = 1.55 ± 0.01 

fcT„ = 1.58 ±0.01 
4 _ « 7 +9 

fcTn = i.43±g:?s 

^4-ray = 12.7_4 


^diskbb = 14.2_ o ; 7 

Aliskbb = 8.3 ± 0.2 

7.71 — 1 1 Q + 016 
"" J ray — l.±0_ q Q3 

Adiskbb = 9.14 ±0.26 
UT — 1 R4+° 14 

K-tray — l.O1_ 01g 

^diskbb = 8.5 ± 3.6 
fcT ray = 1.1 ±0.1 


/dray = 1.68_q 

A/h = 1.46 ± 0.06 
A H =0.8 ±0.1 


17.3/29 
542/500 
1276/1341 
23.2/20 



□Exponentially cut-off power-law of the form Ap^E r exp(— E/Ef a \d) where E is the photon energy and Ef a \^ is the folding energy, both 
measured in keV. 

Epower-law of the form Ap^E~ r , where the symbols have the same meaning as for the cutoffpl model. 

0MCD model after Mitsuda et al. (1984), fcTj n is the temperature at the inner edge of the disk measured in keV (T(r) oc r~ 3 / 4 where r is 
the radius) and the normalization constant is Adiskbb = (fin AO 2 cosi where ri n is the inner disk radius measured in kilometers, d is the source 
distance, measured in units of 10 kpc, and i is the disk inclination. 

0Comptonization model after Titarchuk (1994) for a disk geometry. fcTo is the temperature of the Wien like seed photon input spectrum 
measured in keV, kT e is the plasma temperature, measured in keV, and r is the plasma optical depth. 

Eoptically thin emission spectrum after [Raymond fc Smith (1977)| where fcT ray is the plasma temperature in keV. The normalization constant 
is A ray = 10~ 17 (47rD 2 ) -1 ■ J n c nn&V , integrating over the emitting volume, where n c and njj are the electron and hydrogen particle densities, 
respectively, and D is the distance in cm. 



Note. — The absorbing column, Nn, has been fixed at 1.3 X 10 21 cm -2 for the PCA data, while it was fixed at the value found from the SIS 
data in the case of the GIS data. In all other cases it was a free parameter (measured in 10 21 cm -2 ). The normalization of the SIS1 detector 
with respect to the SIS0 detector was 0.987 ± 0.006, while that of the GIS3 with respect to the GIS2 was found to be 1.043 ± 0.004 in all cases. 
The uncertainties given are at the 90% level for one interesting parameter. 
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Fig. 2. — Same as Fig. [j] for the ASCA data. Left: SIS data. Right: GIS data. Residuals at 1.7 and 2.2 keV are well-known 
ASCA response features. 
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The High Energy Astrophysics Archive (HEASARC) 
contains one ASCA observation of V1408 Aql, made in 
1994 October (Table [l]). A preliminary analysis of these 



data has been presented by Ricci et al. (1996). We ex- 



tracted the data from all four instruments on ASCA, the 
two solid state detectors (SISO and SIS1) and the two GIS 
detectors (GIS2 and GIS3). Due to the uncertainty in the 
intercalibration of the instruments, the GIS and the SIS 
detectors were analyzed separa tely. We describe the data 



extraction process in appendix A. 2 

In Table || we list the results from modeling the data 
with the standard models for the SIS and the GIS, respec- 
tively. The data and residues for the models are shown 
in Fig. ||. Note that due to our extraction procedure 
the model normalizations differ between the detectors. It 
is only pos sible to compare the spectral shapes (see ap- 
pendix |A|). As with the RXTE-PCA data, both the ex- 
ponentially cut-off power-law and the MCD model provide 
a rough description of the data. Due to the higher spectral 
resolution of the ASCA detectors, however, the causes for 
the spectral deviations are now apparent, and include a 
strong deviation at ~ IkeV. We interpret this feature as 
evidence for the presence of line emission at this energy, 
which might come from the iron L complex or emission 
features from other metals suc h as Kg lines from highly 
ionized neon or magnesium (see Nagase et al. 1994). Mod- 
eling the feature with the addition of a simple Gaussian 
line does not result in a markedly improved fit. In contrast, 
including an optically thin thermal plasma spectrum after 
Raymond & Smith (1977) in the spectral modeling process 



results in a dramatic improvement of the fit (A^ cd = 0.65 
for the MCD model) . The best-fit parameters for the ther- 
mal plasma are similar to the disk temperature found with 
the MCD model and the emission line spectrum is domi- 
nated by emission around 1 keV. 

In order to further check whether the 1 keV feature is 
always present in the X-ray spectrum of V1408 Aql we 
turned to archival ROSAT position sensitive proportional 
counter (PSPC) data. The observing log for this observa- 
tion is given in Table [l an d th e data extraction procedure 
is described in appendix [A.3| . As can be seen from our 
fit-results in Table g, the ROSAT data give similar results 
as the ASCA data. In fact, the ROSAT data require the 
presence of the line emission component to provide satis- 
factory fits (see also Fig. ||) ■ 

We note that the PCA data shows weak residuals in the 
region of an Fe line. An MCD model with weak power law 
tail, for example, admits the inclusion of a 6.6 keV line 
with width 0.8 keV and equivalent width 80 eV. Such a 
weak, broad line, however, is comparable to the remaining 
uncertainties in the PCA response matrix, and we there- 
fore cannot be confident of its significance nor of its pa- 
rameters. Adding an Fe line (with energies ranging from 
6.4 to 7.1 keV) to the models of the ASCA data also does 
not significantly improve the fits. Limits to the equivalent 
width of any line in this region were of O(10 eV), which is 
comparable to the equivalent width of the Fe line present 
in the best fit Raymond-Smith models. We note that con- 
trary to the EXOSAT and Ginga data, the ASCA data 
also do not show strong evidence of a hard tail. The up- 
per limit to the flux from a power law tail was 12% of the 
2-10 keV flux in the cutoff power law model of the ASCA 
GIS data. The upper limits for the MCD models and for 



the SIS models were 3-20 times lower. We therefore can- 
not nikj-nrt^thcj^ 7.06 keV line claimed 
by Ricci, Israel & Stella (1995) was associated with the 



presence of a hard tail. 

3. TIMING ANALYSIS 
We employed Fourier techniques, in the same manner 



as for our RXTE observ ations of Cyg X-l ( Nowak et al, 
1999a|) and GX 339-4 (|Nowak et al. 1999b|), to study 



the short timescale variability of V1408 Aql. We use 
the same techniques for estimating dcadti me corrections 
(Zhang et al. 1995; Zhang & Jahoda 1996) to the Power 
Spectral Density (PSD), and for estima ting uncertainties 



and the Poisson n oise levels of the PSD (Leahy et al. 1983 



van dcr Klis 1989) as in our previous RXTE analyses. We 
use lightcurves with 2 -5 s resolution, from the PCA top 
layer data only, over the energy range « 1-7.2 keV (abso- 
lute PCA channels 1-20), in the analysis discussed below. 
We also searched 2 -11 s lightcurves over the same energy 
range for high-frequency features, but none were found 
above the Poisson noise limits. 

As the source intensity did not appear to vary over the 
course of the observation, we created a single PSD. The 
results are presented in Figure || for a normalization where 
integrating ove r positive frequency yields the mean square 
varia bility (sec Belloni & Hasinger 1990; Miyamoto et al 



1992). Note that above / = 10 1 Hz the power is com 



pletely consistent with Poisson noise. We estimate that 
the background contributes 13cps to the lightcurves, com- 
pared to 210 cps for the signal. Based upon these count 
rates and from calculating the PSD of the background 
lightcurve generated using the RXTE software, we find 
that the PSD observed between / = 10~ 3 -10~ 2 Hz is con- 
sistent with background fluctuations. We find that the 
upper limit to the root mean square (rms) variability be- 
tween 10" 3 -16Hz is 4%. 

4. LONG-TERM VARIABILITY 

We used data from the All Sky Monitor (ASM) on 
RXTE to study the long-term behavior of V1408 Aql. The 
ASM is an array of three shadow cameras combined with 
position sensitive proportional counters that provides for 
a quasi-continuous coverage of the sky visibl e from RXTE 
(Levineet al. 1996; Remillard & Levine 1997). Lightcurves 



in three energy bands — 1.3-3.0keV, 3.0-5.0 keV, and 5.0- 
12.2 keV — as well as over the whole ASM ba nd are pub- 
lically available from the ASM data archives (Lochner & 
Remillard 1997| ) . Typically there are several 90 s measure- 
ments available for each day. 

In Figure | we present the ASM data of V1408 Aql that 
were available as of 1998 November 20. The date of our 
pointed RXTE observation is indicated by an arrow in this 
figure. Several features are immediately apparent in these 
data. The count rate light curve shows significant variabil- 
ity with fluctuations up to O(50%) of the mean. These 
fluctuations occur on O(100 day) timescales. The color 
lightcurve (we show the 1.3-3.0keV lightcurve divided by 
the 5.0-12.2 keV lightcurve) shows significantly less vari- 
ability, with peaks in the softness of the source occurring 
on O(400 day) timescales. Furthermore, the peaks in the 
softness of the source seem to be correlated with dips in 
the intensity of V1408 Aql. The features in the light curve 
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appear to be associated with possible long term periodic- 
ities. 

We determined the significance of these possible long 
term peri odicities by computing th e Lomb-Scargle Pcri- 
odogram ( |Lomb 1976| ; |Scargle 1982| ) for the 1.3-12.2keV 
band for 6-day averages of the ASM lightcurves. We aver- 
aged data where the best fit to the source positi on and flux 
in an ASM obs ervation had a xjr od < 1.5 (see Lochner & 
Remi Uard 1997 ) in each of the three ASM energy channels. 
The periodogram presented in Figure || shows evidence of a 
117 day, 235 day, and a 352 day periodicity. Each of these 
periodicities is significant at gre ater than the 95% level, as 



determined by the methods of Home & Baliunas (1986) 



We note that the Lomb-Scargle periodogram does not as 
sume the presence of ha rmonics; this is a result of the anal 
ysis. Epoch f o lding (sec Leahy ct al. 1983 ; |Schwarzcnberg- 
Czern |y 1989| pavics 1990| ) of the ASM lightcurves also 
shows evidence of these periodicities, although each pe- 
riod has uncertainties of approximately ±10 days. The 
evidence for a periodicity in the color lightcurve is some- 
what weaker. Only the longest period appears, with an 
approximately 370 day period, and then only at the 50% 
significance level in a Lomb-Scargle periodogram. 

Figure shows the result of fitting three harmonically 
spaced sinusoids to the count rate and color lightcurves. 
In these fits, the periods were constrained to be within 
a few days of the periods found in the Lomb-Scargle 
perio(|B"gram of the count rate lightcurve; however, the 



phases of the sinusoids were left completely tree. For the 
count rate lightcurve, the amplitudes of the sinusoids are 
0.32 cps, 0.26 cps, and 0.41 cps for the fundamental, first 
harmonic and second harmonic. For the color lightcurve, 
the respective amplitudes are 0.1, 0.05, and 0.01. Further- 
more, the phases of the sinusoids are displaced from those 
of the count rate lightcurves by of 0(ir). In Figure |5| we 
also show the lightcurves folded on the 117 day period. 
Note that the folded color lightcurve indeed exhibits very 
little variation on this timescale. The count rate lightcurve 
shows significantly more periodic structure. The low flux 
points, however, display the most variations from phase 
bin to phase bin. Partly this could be due to interference 
from the 235 and 352 day periods. Additionally, if this 
periodicity is due to inclination effects in a warped disk, 
as we further discuss below, the low flux points come at 
times when the disk is at its greatest inclination to our line 
of sight. The lightcurve is most sensitive at these times to 
small changes in disk thickness and/or shape. 

Long timescale periodicities and quasi-pcriodicities are 



model, and Comptonization. We did not see any evidence 
for a hard power-law tail similar to that seen in previous 
observations. We have also found evidence for a spectral 
feature at ~ 1 keV, which we interpret as emission from the 
iron L complex or as Ka lines from highly ionized metals. 
No short term variability in excess of the noise was de- 
tected from the source, and the upper limit to the rms 
variability between 10 _3 -16Hz is 4%. On long timescales, 
we found evidence for periodic variability on a time-scale 
of about 117 days in the soft X-ray luminosity, and evi- 
dence for a periodic softening of the X-ray spectrum on a 
350-400 day timescale. Below, we discuss interpretations 
of these results. 

SPECTRAL CONSIDERATIONS 

Although the Comptonization model appears to provide 
the best fit to our broad-band RXTE data, we do not 
regard it likely that Comptonization is indeed the physi- 
cal process responsible for producing the X-ray spectrum. 
As we have shown in §||, part of the small Xied obtained 
for Comptonization is attributable to the comparably high 
seed photon temperature, kTo — 0.34 keV, which mimics 
the soft excess seen in the RXTE data. Also, the best fit 
parameters hint at a very cold and optically thick Comp- 
tonizing plasma with an optical depth of almost 10. Com- 
monly assumed models for Compton coronae, such as ad- 
vection domin ated accretion flows (Esin, McClintock & 
Narayan 1997|) or ot her 'sphere plus disk' coronal mod 



relate 'ely common in ASM observations of binary sources 
(RemiUard 1997, private communication). Evidence for a 
294 d periodicity in Cyg X-l has been previously rep orted 
( |Kcmp et al. 1983| ; [Pricdhorsky, Terrell fc Holt 1983| ), and 
is readily apparent in the ASM data during the hard state. 
A 198 day periodicit y also has been obser ved in LMC X-3 
(Cowley et al. 1991; Wilms et al. 1998b), and a possible 
240 day periodicity appears in ASM dat a of the low/hard 
state of GX 339-4 ([Nowak et al. 19991]). 



els (Dove et al. 1997), have considered only hot, optically 
thin to moderately optically thick coronae. It is not clear 
whether a cool and very optically thick corona can be made 
energetically self-consistent, nor is it clear what physical 
processes would lead to such a configuration. We therefore 
conclude that Comptonization is an improbable physical 
mechanism for producing the observed soft spectrum. 

The accretion disk spectrum and the exponentially cut- 
off power-law both provided similar quality fits and yielded 
almost indistinguishable residues. The MCD model, how- 
ever, seems the better phenemonological representation of 
the underlying physical mechanism for producing the ob- 
served spectrum. The best fit parameters for the expo- 
nentially cutoff power-law span a wide range (including 
a negative photon index in the RXTE spectrum), and in 
many ways appear to be "mimicing" the features of the 
MCD model. The MCD models, on the other hand, have 
best-fit spectral parameters that are all similar for each of 
the independent observati ons. More importantly, optical 
and infrared observations ( Cowley, Hutchings fc Crampton 



1988 



5. DISCUSSION — THE NATURE OF V1408 AQL 



_ |Shahbaz et al. 1996| ; [Hakala, Muhli fc Dubus 1999[ ) 
provide independent evidence for the presence of an ex- 
tended accretion disk in V1408 Aql. As we discuss further 
below, additional independent evidence for the assumption 
that the X-rays are dominated by the accretion disk comes 
from the presence of the long term spectral variability. 

We note that the line features apparent in the ASCA 
and ROSAT data are also consistent with an accretion 
disk picture. Line features around 1 keV are a common 
occurrence in photoionized plasmas close to sources emit- 



To summarize, our spectral analysis has provided evi- 



ting h ard X-rays (e.g., in eclipse in Vela X-l, Nagase et al 



dence for a very soft spectrum which can be satisfacto- 
rily modeled with any of the three traditional models used 
here, namely the exponentially cutoff power-law, the MCD 



1994). We would also expect such features in models with 
warped accretion disks similar to those of |Schandl (1996)] 
(see discussion below). Iron L features and Ka lines from 
Mg and Ne are also predicted in models for reflection off 
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ionized accretion disks ( Ross fc Fabian 1993 ), and are in 
fact observed in several NS-LMXB such as Cygnus X-2 



JVrtilek et al. 1986| ; [Kallman, Vrtilek fe Kahn 1989|) , al 



beit the complexity of the observed line shapes makes a 
direct co mparison between the data and the models diffi- 
cult. See Kallman et al. (1996) for a detailed discussion of 
these features. 

LONG TERM VARIABILITY 

The timescales of the periodicities observed with the 
ASM are comparable to the timescales expected from pro- 
cessing accretion disk warps, whether they are dr iven by 



] (jll alju |MaluiiLy, DtgLliiuui Piiuglt 



the radiation [J icssmc iiisLabilily discuii'ii'd jr^PrTrrglT 

(1996 
oney 
1998| 
X-l 



fe DtgLliiMii 199T| , [Maluiit y , DtgLhiiaii fc Nuw<S 



timescales. The warped disk picture could also explain the 
observe d change in the optical lightcurv c recently discov- 
ered by Hakala, Muhli & Dubus (1999) as a precession of 
a warp on long timescales. 

SHORT TERM X-RAY VARIABILITY 

Although we only have upper limits for the amplitude of 
the 10~ 3 -16 Hz variability, these limits are consistent with 
the few observatio ns of BHC and so me NS-LMXB in nearly 
"pure" soft states (Miyamoto 1994). For examples of BHC 
high/soft states, with little or no discernible hard tail 
wherein shor t term variability is presented, s ee |Grebcncv 



(1988 



or by an X-r ay heated wind as for models of Her 



suda 



et al. (1991) for an observation of GX 339—4, Treves et al 



for an observ ation of LMC X-3, Ebisawa, Mit i 
Inoue (1989)| for an observation of LMC X- 1 , anc 



Schandl 1996). As radiation pressure must typi- 



cally strongly dominate gas pressure before a win d can 



be la unched, the former mechanism may dominate (Mai 
oney \z Begelman 1997 ), at least for warps large enough 



such that the outer disk is effectively lluminated by the 
X-ray flux from the inner disk. This radiation pressure 
driven instability is fairly generic, and is expected to cause 
a radiatively efficient (i.e., non-advection dominated) ac- 
cretion disk to warp and precess on 0(100 day) timescales. 
The observed ratio between the precession period and the 
orbital period in V1408 Aql is too long to be exp lained by 
a tidally forced precession of the accretion disk ( Larwooq 
19980 

In a warped disk scenario, the long term modulations 
could be due to a combination of the flux varying as the 
cosine of the inclination angle, as well as due to obscura- 
tion of the inner disk by the outer disk. For the former 
effect, we note that i f the inclination of V1408 Aql is 70° 



75° as suggested by Hakala, Muhli fc Dubus (1999), then 



relatively modest inclination variations of ±10° can yield 
the observed X-ray luminosity variations. The softening 
of the spectrum observed on the 352 day timescale could 
be due to a warp periodically obscuring the central re- 
gions of the accreting system, which would explain why 
we do not detect the hard power-law tail seen by the pre- 

ncb 



VIOUS 



ubsei vat iuiis. In aiialug y Lu uLhei sufL 



sottrcES s 



as L MC X - 3 (Wilms et al. 1008b), we can a33umc that 



this tail is produced in a small and comparably cold ac- 
cretion disk corona close to the compact object which is 
then obscured by the precessing warp. A prediction of this 



scenaiiu, Lheiefuie, is llial a lung Leiin inuiiiluiiiig cain- 
paign ' with an instrument capable of detecting the hard 
power-law tail (e.g., RXTE or BeppoSAX), will detect a 
periodic change in the flux level of the power law, including 
a periodic disappearance of this tail. 

One alternative explanation is that the corona is covered 
by the rim of a (geometrically thick) accretion disk. Un- 
like the warped disk scenario where the relative inclination 
of the disk to our line of sight does not change on orbital 
timescales (see Figure ^), in a disk rim scenario the rim 
is caused by interaction of the accretion stream with the 
outer edge of the disk. Our relative view t hrough the rim 
there fore chang es on orbital timescales (see Hakala, Muhli 
& Duyus 1999 and references therein). This seems to be 
less likely than the warp scenario, however, as contrary to 
the optical and infrared data there appears to be no evi- 
dence for a modulation of the X-ray spectrum on orbital 



[Miyamoto et al. (1994)| for observations of Nova Muscae 
The PSD presented in these works typically have a PSD 
level of (rms) 2 /Hz w 10 -3 at 0.01 Hz, which decreases as 
j -0.7 j? or higri er Fourier frequencies and an rms variability 
of re 3% in the 10~ 2 -30 Hz range. This is slightly below 
the upper limits presented in Figure [I| 

The 'normal branch' of the NS-LMXB GX5+1 also has 
a similar amplitude and shape PSD as described above for 
high/soft state BHC; however, its energy spectrum con- 
sists of both a 1 keV MCD component and a 2 keV black- 
body spectral component (Miyamoto 1994, and references 
therein). The high/soft state o f Cir X-l has been sim- 
ilarly modeled (Miyamoto 1994, and references therein). 
If V1408 Aql had weak 1-10 Hz vari ability comparable to 
that discussed by Miyamoto (1994) for Cir X-l, our ob- 
servations would have detected it. Other soft neutron star 
sources with luminosities of O(10%) Eddington (the ap- 
proximate luminosity of V1408 Aql, if it were a 1.4 M Q 
neutron star given our hypothesis of a highly inclined 
disk), especially the bright atoll sources such as GX13+1, 
GX3+1, GX9+1, and GX9+9, can also exhibit "very low 
frequency noise" with approxi mately 5% rms variability 
( |Hasinger fc van dcr Klis 1989| ). The 10" 4 -100s low am- 
plitude variability in the light curves of these sources has 
been interpreted as intermittent, slow nuclear b u rning on 
the s urface of the neutron star (Bildsten 1993; Bildstcn 
1995). Low frequency variability at such a level is ab- 



sent in V1408 Aql. Furthermore, bright atoll sources often 
show a 0.1-10Hz po wer spectrum in excess of the upper 
limits discussed here ( Hasinger & van dcr Klis 1989|). The 
level of the 0.1-10 Hz PSD seen in GX13+1 ([Homan et aT 



1998), for example, also would have been easily detected 



in the PSD of V1408 Aql, yet was not. 

BLACK HOLE OR NEUTRON STAR? 

The nature of the compact object as of now is not clear. 
The general picture outlined above is similar to that seen 
in neutron s tar X-ray binaries such as Sco X-l, C yg X-2, 
and others. Yaqoob, Ebisawa & Mitsuda (1993) pointed 
out that the normalization of the best-fit MCD model ap- 
pears to indicate that the compact object is a neutron 
star. This argument, however, strongly relies on the as- 
sumed distance to V1408 Aql, for which no compelling 
measurement exists (a lower limit of 2.5 kpc comes from 
the fact that the ASCA and ROSAT measured Ah val- 
ues are consistent with the full galactic column), and also 
relies on the assumption that the accretion disk is seen 
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Line of Sight 




X-ray Heated Surface Orbital Timescales 

Orbit is 9.33 Hours; Warp Precesses on 1 17 Day Timescales 

Fig. 6. — A possible model for the behavior of V1408 Aql (not to scale). Optical modulation is due to partial eclipsing of a warped 
disk. Long term X-ray variations are due to the changing orientation of the warp. The X-ray power law is due to Comptonization 
within a corona. The ~ 1 keV plasma component is due to a wind or possibly due to a disk atmosphere. 



closer to face on. The recent optical and soft X-ray vari- 
ability measurements, however, make a large inclination 
more probable. 

Taking these points into account and assuming for the 
sake of argument a source distance of 7kpc, then the over- 
all flux of V1408 Aql is comparable to that of the high state 
of the black hole candidate GX 339—4, which is a very 
plausible BHC. The upper limits to the high frequency 
varia bility discussed above are consistent with previously 



obser ved BHC power spectra in high/soft states. If transi- 
tions from the hard state to the soft state occur at 5%-10% 



of the Eddington luminosity (see Nowak 1995 ), then the 
compact object in V1408 Aql is consistent with being a 
2-3 M black hole. Thus, although there is no compelling 
evidence that V1408 Aql contains a black hole, there also 
is no compelling evidence that V1408 Aql is a neutron star. 

CONCLUSIONS 

X-ray spectroscopy and the study of both the long term 
and the short term variability of V1408 Aql make a system 
geometry as that depicted in Figure [j] seem likely. A low- 
mass main sequence star serves, via Roche Lobe overflow, 



as the donor for a compact object which is surrounded by 



a large accretion disk which in turn dominates the system 
at all wavelength ranges. The accretion disk is surrounded 
by an optically thin plasma, either in the form of an accre- 
tion disk wind or a stationary accretion disk photosphere, 
which emits the observed X-ray line radiation. A small hot 
corona directly surrounding the compact object produces 
the hard X-ray power-law. The whole accretion disk pre- 
cesses on a time scale of about 117 d, obscuring the central 
region and causing the power-law tail to periodically dis- 
appear and reappear. Also on these long timescales, the 
changing view of the warp causes the orbital modulation 
of the optical light-curve (due to partial obsc uration of the 
outer accretion disk) to vary from sinusoidal ( Thorstciiscr 
1987) to a more complex pattern (Hakala, Muhli & Dubut 
1999[ 

The nature of the compact object in V1408 Aql is still 
ambiguous. We have put forth a hypothesis, however, that 
might explain the observed phenomenology and makes pre- 



dictions that are observationally testable. X-ray monitor- 
ing over the 117 d period with an instrument like RXTE 
or BeppoSAX should reveal whether the X-ray power-law 
tail really does periodically disappear and reappear as pre- 
dicted by our model. Furthermore, if the source is at 10% 
LecM an( i contains a neutron star, then a bout one "Typ e 
I" microburst per day might be expected ( |Bildsten 1995 ). 
This should be easily observable during such a ca mpaign. 
One migh t also hope to find "kilohertz QPO" (vandci 



Klis 1998), as are often associated with atoll sources. For 



these latter two possibilities, however, we note that some 
of the brighter atoll so urces such as GX13 +1 have yet to 
exhibit kilohertz QPO ( Homan ct al. 1998 , and references 
ther ein), and rarely exhi bit Type I bursts (see, for exam- 
ple, Matsuba et al. 1984, and references therein). Finally, 
high spectral resolution observations as will be provided by 
the upcoming new generation of X-ray instruments, such 
as the gratings on the Advanced X-ray Astronomy Facility 
(AXAF) and the X-ray Multiple Mirror Mission (XMM), 
will provide the spectral resolution necessary for resolving 
and studying the Fe L complex. This will allow th e appli- 



cation of p lasma spectroscopic diagnostics (e.g., Licdah] 
ct al. 1992| ) to the study of this fascinating source. 
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APPENDIX 

DATA ANALYSIS METHODOLOGY 
RXTE Data Analysis 

Our RXT E data were analyzed using the same procedure as that for our analysis of the spectrum of GX 339—4 (Wilms 



et al. 1998a). Screening criteria for the selection of good on-source data were that the source elevation was larger than 10°. 
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Data measured within 30 minutes of passages of the South Atlantic Anomaly or during times of high particle background 
(as expressed by the "electron ratio" being greater than 0.1) were ignored. Using these selection criteria, a total exposure 
time of 27ksec was obtained. To increase the signal to noise level of the data, we restricted the analysis to the first anode 
layer of the proportional counter units (PCUs) where most source photons are detected (the particle background is almost 
independent of the anode layer), and we combined the data from all five PCUs. 

To take int o account the calibrat ion uncertainty of the PCA we applied the channel dependent systematic uncertainties 



described by Wilms et al. (1998a) . These uncertainties were determined from a power-law fit to an observation of the 
Crab nebula and pulsar taking into account all anode chains; however, they do also provide a good estimate for the first 
anode layer only since most of the photons are detected in this layer. 

Since V1408 Aql has a comparably small count rate we are able to use the new background model for the PCA that 
was made available by the RXTE Guest Observers Facility (GOF) in 1998 June. The quality of this model was checked 
by looking at high detector channels which are completely background dominated. Although the measured count rate of 
V1408 Aql was at the high end of the applicability of the new background model, the agreement between the model and 
the measured background was good. This is in part due to the fact that V1408 Aql is a very soft source which allows 
greater latitude in using the background model for faint sources. Remaining background residuals were minimized by 
using the XSPEC "corfile" facility which renormalizes the background flux to decrease the best fit Xred- The corrections 
applied to the background flux were on the order of <;1.5%, indicating that at least for this source the background model 
provides a good background estimate. Since the spectrum is completely background dominated above 20 keV, and due to 
the calibration uncertainty below 3keV, we restricted the spectral analysis to the range from 3 to 20keV. 

For the timing analysis, we generated lightcurves from the 'GoodXenon' data. Note that although there are short data 
gaps of 1-4 s duration that are flagged by commensurate jumps in the value of the time coordinate from one data bin to 
the next, there are occasional data gaps where the extraction software generates a continuous series of time bins despite 
the data losses. These data gaps do not appear in lightcurves generated from the 'standard2f data (which is processed 
by a different event analyzer on-board RXTE). These gaps can be recognized, however, in the high time resolution data 
by searching for any sequence, 1 s or greater in length, of time bins with zero count rate. Four such 'unflagged' sequences, 
with 16 s duration each, were found in our data. (Aside from these four 16 s sequences, there were a few instances where 
two 1/32 s time bins in a row would have zero detected counts. The lack of counts in these bins were consistent with 
counting statistics, and we did not consider these to be data gaps.) The power spectra that we presented in Figure || 
were made from continous data segments without internal data gaps. If we include data segments with the unflagged 
data gaps in the calculation of the PSD, we obtain a low amplitude (5% rms) PSD that is flat from 10 _3 -10 _2 Hz and is 
exponentially cutoff at higher Fourier frequencies. In fact, the presence of unflagged data gaps can be deduced from such 
a characteristic PSD shape (Wijnands 1999, priv. comm.). 

ASCA Data Analysis 

We extracted data from the two solid state detectors (SISO, SIS1) and the two gas detectors ( GIS2, GIS3) onboard 
ASCA by using the standard ftools as described in the ASCA Data Reduction Guide (pay et al. 1998j ). The data extraction 
regions were limited by the fact that all the observations were in 1-CCD mode and that the source was placed close to 
the chip edge. To maximize the extraction regions, we chose rectangular regions of 6' x 8' and 6' x 7' for SISO and SIS1, 
respectively. Choosing a rectangular region does not effect the shape of the extracted spectrum; however, the ASCAARF 
ancillary response matrix generator assumes a circular region, so the flux normalization is slightly off (hence the 30% 
normalization differences between the SIS and GIS detectors in Table ||). For the GIS detectors we chose circular regions 
centered on the source each with a radius of 13'. 

The SIS count rate for V1408 Aql is large enough that the central regions of the CCD suffer from pileup (i.e., two or 
more events being registered as a s ingle event). Estimates of the amount of this pileup can be found in the appendix 



presented by Ebisawa et al. (1996). Based upon our measured spectrum and these estimates, we chose to exclude from 
analysis central rectangular regions with dimensions of 4' x 3' and 3' x 3' for SISO and SIS1, respectively. With these 
exclusions, we estimate that pileup will contribute less than 1% of the counts at 10 keV. 



We used the SISCLEAN and GISCLEAN tools (Day et al. 1998), with the default values, to r emove hot and flicke ring 



pixels. As the spectrum of V1408 Aql is very similar to the low flux level of Cir X-l described by Brandt et al. (1996) , we 



filtered the data with the same cleaning criteria outlined in that work; however, we took the slightly larger va lues of 10 



for the minimum elevation angle and 7GeV/c for the rigidity. Also similar to the work of Brandt et al. (1996), we formed 
background estimates by extracting a circular region of radius 5' near the edge of the detector for the GIS observations. 
For the SIS observations, we chose L shaped regions near the corner of the chip opposite from the source. Background, 
however, contributes relatively little to the observations. 

We rebinned the spectral files so that each energy bin contained a minimum of 20 photons. We retained SIS data in the 
0.6 to lOkeV range and GIS data in the 1 to lOkeV range. The cross-calibration uncertainties among the instruments 
were accounted for by introducing a multiplicative constant for each detector in all of our fits. As discussed above, the 
resulting data files showed reasonable agreement between all four detectors. 

ROSAT Data Analysis 

The extraction of the ROSAT spectrum was perfo rmed using the standard R OSAT PSPC data analysis package Ex- 
tended X-ray Scienti fic Analysis System (EXSAS) ( Zimmermann et al. 199g| ) following the procedures described by 



Brunner et al. (1997). Source counts were extracted from a circular region centered on the position of V1408 Aql with a 
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radius of 2', while the background was extracted from an annulus centered on the source from which source counts from 
detected background sources were removed. A correction for the telescope vignetting was applied to the standard ROSAT 
response matrix. The spectrum was then rebinned into 26 channels of ~ 10000 counts each to ensure an even signal to 
noise ratio over the whole ROSAT energy band. As for RXTE and ASCA, the spectral analysis of the extracted data was 
then performed with XSPEC, ignoring data measured below 0.5 keV and above 2.5 keV. 
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